Asymmetric light transmission through a photonic crystal with relaxing Kerr 

nonlinearity 
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The idea of asymmetric light transmission through a photonic crystal with relaxing Kerr non- 
linearity is presented. This idea is based on the symmetry breaking due to the self-trapping of 
intensive pump pulse inside the structure. As a result, the asymmetric transmission of the sec- 
ondary (probe) pulses through the perturbed photonic crystal can be observed. The cases of short 
and long (quasi-continuous) probe pulses are considered. 
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I. INTRODUCTION 

The problem of asymmetric light propagation attracts 
much attention in modern scientific literature. The 
essence of this phenomenon is the situation when light 
transmission through the optical system in two opposite 
directions give different radiation characteristics. In the 
limiting case, one can observe the so-called optical diode 
action when light is transmitted or entirely blocked de- 
pending on the propagation direction. The importance 
of such effects is connected with the perspectives of inte- 
grated optics and all-optical logics and can lead to con- 
struction of all-optical diodes, isolators and other devices. 
To achieve asymmetric or unidirectional transmission, a 
number of optical structures were recently proposed on 
the basis of different operation principles. A number 
of suggestions is connected with the use of multilayer 
systems (or one-dimensional photonic crystals), for ex- 
ample, the passive optical diode based on nonlinear ef- 
fects near the band edge [l], 0] , the hybrid Fabry- Perot- 
resonatorphotonic-crystal structure with strongly asym- 
metric resonant modes [||, photonic crystal heterostruc- 
ture with broken spatial inversion symmetry metal- 
dielectric multilayer diode due to a light-tunneling mech- 
anisni[5j , asymmetric multilayer with left-handed mate- 
rials [6[ and so on. A class of systems utilizes the effects 
of photonic crystal defects to achieve unidirectionality, 
such as coupling of nonlinear defects Ml or microcavi- 
ties with the modes of opposite parity [8(, asymmetri- 
cally placed nonlinear defect in the structure with chi- 
ral @ or magnetic layers [l0| . photonic crystal waveg- 
uide with uncoupled microcavities which is subject to 
a Fano-like mechanism and asymmetric light localiza- 
tion [ll|. In ref. [l2[ two-dimensional photonic crystal 
heterostructure was used to obtain diode action due to 
a self-collimation behavior. Another type of structures 
considered in literature uses the so-called metamateri- 
als, for example, planar split-ring resonators with chi- 
ral symmetry breaking [13j . bilayered chiral metamate- 
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rial with polarization-dependent asymmetry |14 |. artifi- 
cial structure with polarization selectivity due to magne- 
toelectric coupling and light tunneling [15]]. Asymmetry 
in transmission was recorded in the nonlinear waveguide 
with a localized phase discontinuity 16], nonlinear di- 
rectional couplers with longitudinally varying coupling 
coefficient [17j], etc. A similar effect of asymmetric re- 
flection suppression (known as unidirectional invisibil- 
ity) was found in parity-time symmetric periodic struc- 
tures [H GJ]. Finally, it is worth to mention the pos- 
sibility of asymmetric transmission of ultrashort pulses 
through a resonantly absorbing two-level medium dis- 
cussed in ref. (20j . 

In this paper we theoretically search for evidence of 
asymmetric transmission of ultrashort light pulses inter- 
acting with a one-dimensional photonic crystal (multi- 
layer structure) made from materials possessing relaxing 
cubic (Kerr) nonlinearity. The need to take noninstanta- 
neousness of nonlinearity into account was realized soon 
after appearance of nonlinear optics as a separate disci- 
pline. To name a few examples, one can mention studies 
of relaxation effects on laser beam self-focusing [2l| or 
parametric amplification [22j . The effects of finite relax- 
ation time become apparent when this time is comparable 
with pulse duration. This situation seems to be very ac- 
tual nowadays as the methods of production of ultrashort 
(femtosecond) pulses are established. In the context of 
photonic crystals, the concept of relaxing Kerr nonlin- 
earity was applied to the phenomena such as pulse com- 
pression [23| and optical switching (24|. The latter work 
also suggests that the polymeric materials are perspec- 
tive media for applications due to both high nonlinearity 
coefficient and fast relaxation. 

The main feature of the nonlinear photonic structure 
important for our study is the self-trapping effect ana- 
lyzed in our previous works [25l - l27j . This effect occurs 
when both periodic modulation of linear refractive index 
and relaxation of nonlinearity are present. In this letter 
we suggest that self-trapping can be used to break the 
symmetry of the photonic crystal: the powerful pump 
pulse trapped inside the system changes the distribution 
of the refractive index which, in turn, influences the inter- 
action of the secondary (probe) pulse with the structure. 
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Such a variant of pump-probe technique is possible due 
to the relaxation of nonlinearity and allows to use pump 
and probe one after another. Simultaneously, this solves 
the common problem of separation between pump and 
probe signals. The paper includes the formulation of our 
theoretical model, the description of the main features 
of self-trapping effect, the results on asymmetric trans- 
mission in the cases of short-pulse and quasi-continuous 
regimes. Finally, we discuss the importance of probe in- 
tensity to obtain the asymmetric transmission. 



II. THE DESCRIPTION OF THE MODEL 

In this letter we consider an ultrashort (femtosecond) 
pulse interaction with nonlinear photonic crystal taking 
into account the relaxation process. Light propagation in 
this system is governed by the Maxwell wave equation, 



d 2 E 1 d 2 {n 2 E) 



dz 2 



dt 2 



0. 



(1) 



where E is the electric field strength, c is the light speed, 
n = n(z,t) is the nonlinear refractive index. The relax- 
ation of nonlinearity is described by the Debye model, 



Til 



dS n(I,t) 
~~dt 



+ 5n(I, t) = n 2 I, 



(2) 



where Sn is the nonlinear part of refractive index, n 2 
is the Kerr nonlinear coefficient, / = \E\ 2 is the light 
intensity. The full refractive index is 



n(z, t) — n (z) + Sn(I, t), 



(3) 



where n (z) is the linear part of the refractive index vary- 
ing periodically along the z axis and defining the struc- 
ture of the photonic crystal. To simulate pulse propa- 
gation in this system, we solve self-consistently together 
with the eqs. ((2]) and ([3|) using the finite-difference time- 
domain method described in ref. 25] . In our calculations 
we assume the light pulses to have the Gaussian shape 
with the amplitude A = A m exp(—t 2 /2t 2 ), the duration 
t p = 30 fs, and the central wavelength of the initial pulse 
spectrum A c = 1.064 /im. As regards the photonic crys- 
tal (see the upper panel of fig. [J), the linear parts of 
refractive indices of the alternating layers are n a = 2 
and rib = 1-5, respectively; their thicknesses a — 0.4 and 
b = 0.24 /im; the number of layers N = 200. The nonlin- 
ear coefficient of the materials is defined through the non- 
linear term of the refractive index, so that 712/0 = 0.005; 
this means that the pulse amplitude is normalized by the 
value Aq = t/Iq. The relaxation time of the nonlinearity 
of all layers is t n i = 10 fs (fast electronic Kerr nonlinear- 
ity). Note that wc adopt equal nonlinearity parameters 
for both types of the layers that do not limit generality 
as was shown in the study of the different variants of the 
structure including the half-linear one [27| . The changes 
in light spectra in such situation were studied in ref. [26j , 
so that we do not consider the spectral effects in this 
paper. 
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FIG. 1. (Upper panel) The scheme of the one-dimensional 
photonic crystal. (Lower panel) Dependence of the output 
light energy (normalized to the input energy) on the peak 
amplitude of the incident pulse. The energy is integrated 
over the time 300£„. 



III. SELF-TRAPPING EFFECT AND TRAP 
POSITION 



The essence of self-trapping inside the photonic crys- 
tal is the localization of light pulse due to formation of 
a nonlinear "cavity" (or trap) which dynamically adjusts 
its properties to light intensity and which is stabilized 
by nonlinearity relaxation. This trap has a form of bell- 
shaped distribution of intensity and refractive index vari- 
ation. This variation preserves its form for a relatively 
long time (more than 100 transmission times in the lin- 
ear case) [23]. Thus, there are three temporal parameters 
of the system: pulse duration t p , nonlinearity relaxation 
time t n i, and pulse propagation time t prop . Obviously, 
the effects of relaxation become important when t p and 
t n i are comparable; t prop must be not too short for relax- 
ation processes to have enough time to show themselves. 
In other words, the structure must be long enough for a 
pulse to be trapped inside it. For our photonic crystal 
described above, t prop ~ 30i p which is a large enough 
value. 

The effect of self-trapping can be indicated by the value 
of light energy leaving the photonic crystal as a result 
of reflection and transmission. This value is calculated 
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FIG. 2. Distribution of light intensity inside the photonic 
crystal (at t — 200t p ) for different pump pulse amplitudes. 
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FIG. 3. Intensity profiles of 10 probe pulses transmitted from 
left to right (LR) or vice versa (RL) in absence of the pump 
pulse. The amplitude of the probe pulses is A' m — 2Ao, the 
interval between the pulses is 10t p ; the first pulse starts at 
t = 200t D . 



by the integration of radiation intensity over the certain 
time interval. If this output energy normalized by the 
energy of the initial pulse and integrated over a large 
enough time is much less than the unity, one can claim 
that the pulse is trapped inside the structure. The re- 
sults of calculation of the output energy as a function of 
the incident pulse amplitude for the above stated param- 
eters are shown in fig. [T] This value is integrated over 
the time 300i p which is about 10 times larger than the 
propagation time of the pulse through the photonic crys- 
tal in linear regime. It is seen that for the amplitudes 
larger than the critical one (about 2.2Ao), the output 
energy abruptly drops which is the evidence of the pulse 
self-trapping. There is the region of optimal light-matter 
coupling. However, as the amplitude increases, the out- 
put energy goes up due to reflection. This is explained 
by the trapping of the pulse near the very entrance of the 
structure [25j . 

As it was mentioned above, the result of the pulse self- 
trapping is the formation of the quasi-stationary trap 
confining most part of the pulse energy. This trap can be 
viewed as a characteristic distribution of light intensity 
inside the system (its total length is 128 /im). In fig. [2] 
we plot such distributions at the instant t = 200i p after 
the pulse incidence and at different values of the pulse 
amplitude. It is seen that at low amplitudes (A rn = Aq 
and 2Aq), the intensity of radiation in the photonic crys- 
tal is very low and its distribution looks disordered and 
is extended over the whole structure. Above the critical 
amplitude (A m = 3Aq and so on), the distribution takes 
a bell-shaped form (due to the similar shape of the pulse) 
which is formed closer to the entrance of the system for 
stronger pulse amplitudes. This change of position is 
connected with the increase of reflection with the am- 
plitude noted in fig. [1] Simultaneously, the peak value 
of intensity inside this trap tends to decrease (the trap 



gets wider) though its absolute value can sometimes be- 
come larger due to a larger absolute value of the pulse 
amplitude. 

The most important feature of the distributions con- 
sidered is their position inside the photonic crystal. The 
initial pulse propagates from left to right, so that the po- 
sition of the trap is closer or farther from the left side of 
the structure depending on the pulse amplitude. Obvi- 
ously, the photonic crystal which was initially perfectly 
symmetric obtains some attribute of asymmetry when 
the trap is located inside it. Let us calculate how this 
asymmetry shows itself in the secondary (probe) pulse 
propagation in two variants: from left to right (LR) and 
vice versa (RL). The initial pulse forming the trap is fur- 
ther called the pump pulse and is launched at the instant 
t = (the time of maximum incidence at the input). 



IV. SHORT PROBE PULSE TRANSMISSION 

To study the problem of asymmetry, we take a series of 
10 probe pulses: the first starts at the instant t — 200i p , 
while the interval between the other is 10t p . All probe 
pulses are identical and have the same subcritical ampli- 
tude A' m = 2Aq. First of all, we consider propagation 
of LR and RL pulses in the absence of any pump pulse, 
i.e. without a trap inside the structure. Figure |3] shows 
the resulting profiles of 10 subcritical probe pulses trans- 
mitted through the nonlinear photonic crystal. One can 
easily see that the transmission in this case is perfectly 
symmetric: the profiles are identical for both directions, 
though every individual pulse behaves in a different way 
comparing with its neighbors. This is the evidence of 
interaction between the pulses which occurs due to the 
shortage of the interpulse interval. 
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FIG. 4. (a) Intensity profiles of 10 probe pulses transmitted 
from left to right (LR) or vice versa (RL) for different am- 
plitudes A m of the pump pulse forming the trap, (b) The 
relation of the peak intensities Ilr/Irl of the 1st, 4th, 7th, 
and 10th probe pulses versus the amplitude of the pump. The 
parameters of the probe pulses are the same as in fig. [3] 



Now we can analyze the influence of light distribu- 
tions formed by the initial (pump) pulse (see fig. [2]) on 
the properties of the transmitted probe pulses. The re- 
sults of calculations of intensity profiles at different am- 
plitudes of the pump pulse are depicted in fig. |U(a). The 
probe pulses start at t = 200i p as previously. It is seen 
that even at the subcritical amplitude of the pump pulse 
(A m = 2Aq), when the residual light intensity inside the 
photonic crystal is very low (I ~ 10~ 6 io)i the perfect 
symmetry between LR and RL transmission is already 
broken. Of course, the difference between these two cases 
is not dramatic, though it gets more pronounced for the 
latest probe pulses than for the first ones. This is the case 
for all the variants shown in fig. |U(a) and is perhaps con- 
nected with the influence of the first pulses on the trap. 
However, the correspondence between the LR and RL 
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FIG. 5. Intensity profiles the quasi-continuous probe pulse 
transmitted from left to right (LR) or vice versa (RL) in ab- 
sence of the pump pulse, (a) The case of the linear photonic 
crystal; (b) nonlinearity without relaxation; (c) relaxing non- 
linearity. The amplitude of the probe pulse is A' m — 2Ao, its 
duration is 1 ps; the pulse starts at t = 200t p (as previously, 
t p = 30 fs). 



pulses is still close for A m = 4Aq. For larger amplitudes 
(A m = 7Aq and 9Aq), as the trap forms closer to the en- 
trance of the system, the difference appears even between 
the first probe pulses moving in opposite directions. They 
are not equal not only in intensity, but also they need sig- 
nificantly different times to transmit through the struc- 
ture: LR pulses propagate faster than RL ones. The 
disparity in intensity grows for latest pulses, so that they 
can differ by a factor of two and more. This point is 
illustrated by fig. @|b) where the ratio of the peak in- 
tensities Ilr/Irl is shown. It is seen that the deviation 
of this value from unity tends to get greater for the lat- 
est probe pulses and for larger pump amplitudes. This 
figure also suggests the choice of proper pump intensity 
to achieve the optimal contrast between the LR and RL 
transmission. 



V. QUASI-CONTINUOUS REGIME 

In the previous section, we have considered the case of 
probe pulses of the same duration as the pump (trapped) 
pulse. Such situation when the probe pulses effectively 
experience the relaxation of nonlinearity (due to t p ~ 
t n i) should be discriminated from the quasi-continuous 
regime. In this last case t p t n i, so that the influ- 
ence of the relaxation on the pulse propagation should 
be not so dramatic. Let us analyze propagation of the 
probe pulse with the duration t' = 1 ps, while leaving 
the relaxation time (t n i = 10 fs) and pump pulse duration 
(t p = 30 fs) unchanged. The results of calculations in ab- 
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FIG. 6. Intensity profiles of the quasi-continuous probe pulse 
transmitted from left to right (LR) or vice versa (RL) for 
different amplitudes A m of the pump pulse forming the trap. 
The parameters of the probe pulses are the same as in fig. [5] 
The duration of the pump pulse is still t p = 30 fs. 
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FIG. 7. The value \Ilr/Irl — 1| for the quasi-continuous 
probe pulse transmitted from left to right (LR) or vice versa 
(RL) versus the amplitude of the probe. The insets show the 
probe pulse profiles for the cases of relatively large and low 
probe intensities. The pump pulse has amplitude A m — 3Ao 
and duration t p = 30 fs. The other parameters are the same 
as in fig. [5] 



sence of pump are plotted in fig. [5] for different cases: the 
linear photonic crystal, the nonlinear one without relax- 
ation, and, finally, the nonlinear structure with relaxing 
Kerr nonlinearity. As previously, we see the symmetry in 
pulse transmission in all cases. It is also worth to point 
to the narrow peaks appearing in the profiles of the pulse 
transmitted through the system with nonlinearity. Obvi- 
ously, this is the clear evidence of the well-known effect 
of modulation instability which leads to the breakup of 
the long pulse into a train of short ones [28| . As the com- 
parison between figs.JSJb) and 03c) shows, the relaxation 
of nonlinearity promotes the modulation instability and 
precipitates the breakup. The study of these processes is 
beyond the scope of this paper, therefore we only make 
here this short comment and turn to discussion of the 
influence of the pump pulse on the probe behavior. 

The dramatically increased efficiency of this influence 
can be easily seen in fig. El Even a subcritical pump pulse 
{A m = 2Aq) appears to be enough to break the symmetry 
in the propagation of the quasi-continuous probe pulse. 
Perhaps, this is connected with the heightened sensitivity 
of such nonlinear processes as the modulation instability 
to the changes in external conditions. This results in the 
strong change of the pulse profile, which, together with 
the asymmetric distribution of the refractive index, leads 
to the strengthening of the dependence of transmission on 
the propagation direction. As a result, the sharp rise-up 
portion of the LR pulse is obtained with the subsequent 
gradual decrease in intensity. The profile of the RL pulse 
is fundamentally different. However, at larger amplitudes 
of the pump pulse (A m = 7Aq and 9^4o)> the situation is 
inverted: the sharp peaks for the RL pulse and unstruc- 
tured radiation for the LR one. This can be explained as 
a result of the probe pulse interaction with the trap local- 



ized near the very left end of the photonic crystal: under 
these conditions, the LR probe is perturbed in the very 
beginning of propagation and does not have enough time 
to form a characteristic profile. The contrast in light in- 
tensity also increases at higher amplitudes of the pump, 
so that the peak intensity in one direction exceeds the 
same value in the other by the factor of three as can be 
seen in the lower panels of fig. 



VI. IS THE INTENSITY OF PROBE 
IMPORTANT? 

This question should by studied to understand, 
whether the self-trapping effect is the only process re- 
sponsible for asymmetric propagation or the interaction 
of the probe pulse with nonlinear structure must be 
taken into account as well. Previously we have consid- 
ered probe pulses with relatively strong amplitude (2Aq), 
not far from the amplitude of self-trapping. Since we 
are interested only in absolute value of asymmetry, we 
calculate the quantity \Ilr/Irl — 1| as a characteristic 
of direction dependence. The results of calculations for 
the quasi-continuous f-ps-long probe pulse are shown in 
fig. [3 One can see that the value of asymmetry strongly 
depends on the probe amplitude and is close to zero al- 
ready for A' m = 0.5Aq. The insets visually demonstrate 
the difference in pulse profiles in the cases of strongly 
asymmetrical and almost symmetrical propagation. Sim- 
ilar dependence takes place for short probe pulses (not 
shown here). 

The results of fig. [7] imply that the strong interaction 
of the probe pulse with the nonlinear photonic structure 
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is necessary to observe asymmetric transmission as well 
as additional refractive index modulation due to self- 
trapping. In other words, strong amplitudes of both 
pump and probe are needed to obtain the described ef- 
fect. It is hard to demarcate the influence of pump self- 
trapping and probe interaction in the situations consid- 
ered. Perhaps, the study of change of asymmetry as a 
function of slow temporal change of the trap can give 
some additional clues. However, the situation seems to 
be similar to that described in ref. [3| where the asym- 
metric distribution of refractive index is used to localize 
the radiation which further nonlinearly interacts with the 
structure. In our case, the asymmetry in refractive in- 
dex is caused by the pump pulse (quasi-stationary trap), 
then the probe pulse nonlinearly interacts with the mod- 
ified periodic system. We should also emphasize that the 
nonlinearity relaxation is necessary to obtain the trap 
[25j and, hence, to obtain a quasi-stationary (and large) 
change of materials refractive indices which influence the 
probe pulse. It is also worth noting that high probe in- 
tensity makes the proposed scenario of asymmetric trans- 
mission different from a typical pump-probe technique. 



trapping in the photonic crystal with relaxing Kerr non- 
linearity. It is shown that the position of the trap formed 
by the high-intensity pump pulse strongly influences the 
propagation of probe pulses and leads to the dependence 
of the transmission characteristics on the propagation di- 
rection. In the case of a series of short probe pulses, this 
asymmetry is sharp enough only for the late components 
of the series, while the single long (quasi-continuous) 
pulse feels the pump radiation much stronger. The light 
dynamics in this latter case are more complex and include 
the processes such as the modulation instability. As a re- 
sult, the efficiency of asymmetric transmission is higher 
under the quasi-continuous conditions than in the short 
pulse regime. The dependence of this effect on probe 
pulse intensity allows to suggest that the role of pump 
comes to formation of asymmetric refractive index modu- 
lation which modifies nonlinear interaction of probe with 
the material of photonic crystal. Additional evidence for 
the asymmetric transmission is likely to be given in the 
study of the effects such as asymmetric bistable response. 



VII. CONCLUSION 

To sum it up, we have considered the possibility of 
asymmetric light transmission in the regime of pulse self- 
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